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Abstract 

We consider the problem of the driven harmonic oscillator in the probability representation of quantum 
mechanics, where the oscillator states are described by fair nonnegative probability distributions of po- 
sition measured in rotated and squeezed reference frames in the systems phase space. For some specific 
oscillator states like coherent states and nth excited states, the tomographic-probability distributions 
(called the state tomograms) are found in an explicit form. The evolution equation for the tomograms 
is discussed for the classical and quantum driven oscillators, and the tomographic propagator for this 
equation is studied. 
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1 Introduction 



Recently [l}|4| the new formulation of quantum mechanics based on identification of any quantum 
state with the tomographic-probability distribution containing the same information on the state that 
is given by the state wave function [5] or the density operator |6||7|was introduced. This formulation 
is the result of many previous attempts to find such a notion of the quantum state which is similar to 
the notion of classical state associated with the probability density of the classical system in its phase 
space (q,p). These attempts provided the notion of the Wigner function W(q,p) introduced in (8), the 
Husimi function Q(q,p) introduced in (91, as well as Glauber function P(q,p) uO\ and Sudarshan function 
ip(z) |11|. All these functions were associated with system density operator by invertable maps. Also 
all these functions are similar in some aspects to classical probability distribution f(q,p)- Nevertheless 
all these functions, called quasiprobability distributions or quasiprobabilities, are not fair probability 
distributions in the phase space. 

At the same time, optical tomogramm w(X, 8), introduced as a technical tool 12fl3 and related to the 



Wigner function by the integral Radon transform [T3] is a measurable in quantum optics experiments [15] 
probability distribution of the photon homodyne quadrature X, where the angle 9 is the local oscillator 
phase. Optical tomogramms were used to reconstruct the Wigner function (considered as a concept of 
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quantum state) from the experimental data. In [T] the notion of the symplectic tomogramm w(X,p,u) 
was suggested. It is also a probability distribution of the photon quadrature X, related with optical 
tomogramm by an invertible map. Also it became clear that tomographic probabilities can be considered 
as primary objects (not as technical tools) and can be identified with quantum states. These tomogramms 
are alternative to the wave function and density operator since they contain complete information about 
the quantum states. 

The problem of driven oscillator is an interesting problem studied first in [l6 - 18 . The coherent 



states of the oscillator and dynamical invariants [19] of such system were found in 20-22 . On the other 
hand, a detail consideration of this problem in the tomographic probability representation has never been 
performed. The aim of our work is to study driven harmonic oscillator in the tomographic probability 
representation. We obtained the symplectic and optical tomogramms as solutions of the tomographic 
evolution equation for coherent and excited states of driven oscillator. 
This paper is organized as follows. 

We find the Green function for the Schrodinger equation in Sec. 2 and the evolution equations for 
the classical probability distribution in Sec. 3, where we also obtain the Wigner function and classical 
and quantum tomograms. We study symplectic and optical tomograms of coherent and excited states in 
Sees. 4 and 5, and finally we give our conclusions in Sec. 6. 

2 Green Function of Driven Harmonic Oscillator 

We study the problem of a quantum oscillator with a time-dependent external force. The Hamiltonian 
of such a system reads (we put m = = 1, H = 1): 

#=f + f-/(*)$ (1) 
Our goal is to find the expression for Green's function G(x,x',t) of Schrodinger evolution equation. For 



this, we employ the method of integrals of motion 19 -22] We consider a classical analog of the Hamilto- 
nian ([l]): 

H(q,p,t) = £ + £-f(t)q (2) 
The classical equations of motion determining the oscillator position and momentum read 

q=p (3) 
P=~q + f(t) (4) 

and the classical trajectories in the oscillator phase space under the initial conditions q(0) = qo,p(0) = Po 
are given by formulae: 

Q{t) = Qo cos t + po sint + / sin(t — r)/(r)rfT (5) 

Jo 

p(t) = po cos t — qo sint + / cos(i — r)/(r)dr (6) 

Jo 

Now we consider the system of equations ^ and ^ as an algebraic system for unknown initial position 
qo and momentum po, respectively. The variables q,p, t are taken as parameters. The solution of the 
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system is given as: 

qo(q,p,t) = qcost — psint + / sin t/(t)c?t = qcost — psini + F(t) (7) 

Jo 

Po(q,p,t) = pcost + qsint — / cos r/(r)dr = pcost + qsint — J(t) (8) 

Jo 

We define operators acting in the Hilbert space of the quantum driven oscillator states as follows: 

% (q, p,t) = q cos t — p sin t + F(t) 
po (q,p,t) = q sin t + p cos t — J(t) 

Calculating the total derivative of the operator qo(q,p,t) with respect to time t we obtain: 

dqo dq 



(9) 
(10) 



dt 



dt 



+ i 



H,qo 



F(t) - f(t) sint = 



(11) 



where we take into account the definition of function F(t) = Jq sinr/(r)dT given by the ([7]). 
Similarly, the total time derivative of the operator po(q,p, t): 

dpo dpo 



dt 



dt 



H,po 



-J(t) + f(t) cost = 



(12) 



where the definition of function J(t) = j cos r/(r)dr from mh is taken into account 



Thus, operators ^ and ( 10 ) are integrals of motion and correspond to initial position and momentum. 



Then these operators must satisfy the equations for the Green's function G(x,x',t) 

%(x)G(x, x', t) = q( x /)G(x, x', t) 
Po(x)G(x,x',t) = -p( x ,)G(x,x',t) 



20 



22 : 



(13) 
(14) 



where the operators in the left-hand sides of equations act on variable x, and in the right-hand sides on 
x' . Now we write Eqs. (13) and (14) explicitly: 



xcosi + i^- sint + F(t) ) G(.r. :v' . I ) — .r'CH.r. x' . I ) 



( d \ 9 

x sin t — i-— cos t — Jit) G(x, x , t) = i——G(x, x , t) 
\ ox J ox' 

where we still need to define function x(t) to obtain more convenient equations, 

F(t) cost - J(t) s'mt = - [ sin(t - r)/(r)dr = -x{t) 

Jo 

It is worth noting that x{t) is solution of the classical equations of motion pi), Q with zero initial 
conditions: q(0) = 0,p(0) = 0. 

Finally, system of equations for defining Green's function G(x,x',t) is: 



(15) 
(16) 

(17) 



dG 
dx 
dG 
dx' 



xcost — x' F(t) 



i sin t i sin t 
x' cost — x x(t) 



i sint 



zsint 



G 
G 



(18) 
(19) 
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where we skipped the arguments of the Green function. 



Now one can integrate (18) with respect to variable x and (19) with respect to variable x' . Afterper- 



forming the integration, we arrive at 



G(x, x , t) = C(t) exp 



+ x' 2 ) cos t - 2xx' xF(t) + x'x(t) % 



2i sin t 



isini 



(20) 



Finally, following the procedure adopted from [19-22 , we obtain the equation for the function C(t) as 
follows: 

,dG (x,x',t) ~ , . 

i -r^- 2 - = H (x) G{x, x', t) (21) 

After some algebra, we obtain an equation that does not contain the variables x, x 1 is obtained: 



ctgt F 2 (t) 



2i 2 sin 2 t 



(22) 



Equation (22) can be simply integrated with respect to time t and one obtains the function C(t) in a 
form containig a constant C: 



C 

= ~i= ex P \ ■> 



1 F 2 (t) 



-At 



o sin r 



(23) 



To calculate the constant C we require that the following condition at the initial moment of time t = 
be satisfied: 

G(x,x',t = 0) = S(x-x') (24) 

After calculating C using ( 24 ) and substituting it into ( 23 ) and ( |20[ ) we arrive at the final expression for 
the Green function of the Schrodinger evolution equation for the wave function, 



G(x, x , t) 



1 



sin t 



exp 



1 F 2 (t) 
sin 2 r 



dr x 



x exp 



+ x' 2 ) cos t - 2xx' xF(t) + x'x(t) 



2i sin t 



is'mt 



(25) 



We consider some particular cases as examples. 

We start with the case of a constant external force, /(t) = fo = const. In this particular case 
we calculate functions x(t) and F(t) according to their definitions (17) and Q correspondingly and 
substitute the result into the final expression (|25|). After some algebra, we obtain 



G{x,x',t) 



1 



\/2m smt 



exp 



\ f 2 t ((x - f ) 2 + (x' - fof) cost - 2 (x - fo) (x' - fo) 



2i sin t 



(26) 



Taking into account that Hamiltonian in this particular case differs from Hamiltonian of a harmonic 
oscillator with shifted equilibrium just in a constant term: 



H=l(p 2 + (q-foY 



fo 
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we write the following set of equalities: 



Hr x \G(x, x' , t) = exp 



.0. 

dt 



G(x, x , t) exp 



f 2 8 
J -£-G(x,x',t) = i—G{x,x',t) 



This means that Green function satisfies the Schrodinger evolution equation. 
Now we considerthe limit case to — > 0. 

First of all we should recover m,uj,h in formulae according to the following rules: 



t — > tot; x 



muj f 

— x;f -> — -=== 
n tovmhoj 



After all the substitutions we have the Green function in the form 
G(x, x , t) = 



m 



tO 



2mh V sino;t 



exp 



2uj 3 mh 



fo"t + 2/ 2 



COS (jjt—1 



sinwt 



2fomu [x + x') 



A cos tot — 1 



sin bjt 



2, ,4 1 



+ m us 



+ x /2 ) coswt — 2xx' 
sin cot 



At 



the Green function reads 
G(x,x',t) = 



in 



2niht 



exp 



M /o * /o* ( x + x ') m(x— x') 



/\2 



m 24 



+ 



+ 



2/ 



(27) 



(28) 



By direct substitution into the Schrodinger evolution equationwe can check the correctness of our cal- 
culations. It is also worth noting that the Green-function phase is the classical action S(x,x',t), which 

i 

satisfies the Hamilton- Jacobi equation, multiplied by -. 

n 

3 Time Evolution of Classical and Quantum Distribution Functions 

Now we consider the propagators of different distribution functions of the system with Hamiltonian 
in the quantum case and Hamiltonian Qin the classical case. We start from the classical distribution 
function and the Liouville equation, 



^+p| + (/(*)-9)|]p( 9 ,P,t) = 



(29) 



Noting that the system trajectories in the phase space are given by Eqs. ([5]), ^ (and, correspondingly, 
integrals of motion are qo(q,p,t),po(q,p,t) from ([7]), Q), we conclude that the evolution of the function 
p(q,p, t) reads: 

Piq^P-, t) = Po(qcost — psint + F(t), qs'mt + pcost — J(t)) (30) 

where po(q,p) is the distribution function at initial moment of time t = and it depends on two variables 
q,p only. 

Now we consider the Moyal's equation [23] for the Wigner's function W(q,p,t): 











dt ^^dq 



U[q^q + 



i d 
2dp 



c.c. 



W(q,p,t) = 



(31) 



5 



where the quantities in the square brackets are pure imaginary and in the equations, pure real ones. The 
term in the square brackets can be rewritten, in view of the fact that the potential energy is given by 
x 2 
~2 



U(x) 



f(t)x, as follows: 



U[q + 



i d 
2dp~ 



(32) 



We are interested just in imaginary part of this expression. After substituting this into (31) we obtain 
the equation: 



dt 



(33) 



It is obvious that Eq. (33), being the evolution equation of the Wigner function W(q,p, t) fully coincides 
with the equation (29) of evolution of the classical distribution function. The dependence of the Wigner 

(34) 



function on time is given by an expression similar to (|30j) , 

W(q,p, t) = Wo(q cost — psint + F(t), qsint + pcos t — J(t)) 



Now we move to the tomographic-probability-distribution functions. 

We start from classical tomogramm distribution. In thisl case we have [4,24 the following evolution 
equation: 




w(X, ji, v, t) = 



(35) 



where F(q, t) has the meaning of acting force and in our case it is given by the expression F(q, t) 
—q + f(t). So, evolution equation in this case has following form: 

dt 



d d ,. . d 



di 



dX 



w(X, n, v, t) = 



(36) 



The integrals of the equation (36) can be easily obtained. They read: 

v o (X, fjL,v,t) = v cos t + fi sin t 
fio(X,fj,, i>,t) = //cost — z^sint 

Xq(X, (J,,U,t) = X — /jLx(t) — vp{t) 



(37) 
(38) 
(39) 



where functions x(t),p(t) are solutions of classical equations of motion ([3]), Q with zero initial conditions. 
One should note that at the initial moment of time t = 0: 



v = v,fJ, = fi,X = X 
hence, evolution of the tomographic-probability distribution is given by: 

w(X, fi, t) = wq {X — fj,x(t) — vp(t), fi cos t — v sint, v cost + [Asint) 



(40) 



where wo(X, fj,, v) is the tomographic-probability distribution at the initial moment of time and it depends 
just on three variables X, fi, v. 

We move now to the quantum tomogram and write the evolution equation according to [1,25 : 

-l 



d_ 
dv 



U\q 



d_ 
dX 




c.c. 



w(X, ix, u,t) = 



(41) 
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d 

We are interested just in the imaginary part of the expression in the square brackets. It reads v— — h 
d 

isf(t)——, and, finally, the evolution equation is 



|-4 + ,| + ,/(t)^)Mx, M ,,,t) = o 



(42) 



that fully coincides with the equation in the classical case (36). Therefore, evolution of quantum tomo- 
gramm is given by the same equation (40). 



4 Symplectic Tomograms of Coherent and Excited States 



Now we study two examples of evolution of quantum tomogramm distributions. 

The first example is the coherent-state evolution, and the second example is the evolution of the nth 
excited state of the harmonic oscillator. 

The wave function of the harmonic-oscillator coherent state [10 11,26] at the initial moment of time 
can be written as follows: 



1 



: exp 



(x - XqY 



+ ipox 



(43) 



We calculate the tomographic-probability distribution of the harmonic-oscillator coherent state at the 
initial moment of time using the following formula (24]: 



w(X,{i,u) 



1 



2vrM 



< '(.</) <-*P ( - % — // { ''hi 



(44) 



After performing the integration, we arrive at 

w (X,fi,u) = 



exp 



(X - vpo - hxqY 

V 2 + /X 2 



(45) 



Now the time dependence of the tomographic-probability distribution is given, in view of (40), as follows:: 

\2\ 



w(X,n,v,t) 



1 



7Tv/ (I 2 + V 2 



exp 



(X - HXd{t) - vpd{t)Y 



V 2 + fj 2 



(46) 



where the functions x c i(t),p c i(t) are solutions of classical equations of motion ([3]), Q under the initial 
conditions x(0) = xo,p(Q) = pq. 

Now we consider the evolution of the nth excited state of the harmonic oscillator with the wave 
function 



^n(x) 



1 



(47) 



2 d n 2 

where H n (x) = (— l) n e x ~^^ e ~ x ^ s Hermite polynomial. We show two ways of finding the evolution of 
the system — using the propagator for the wave function (Green function G(x,x',t) from (25)) and the 
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evolution of the quantum tomogram given by ( 40 ) . 

To study the evolution of the wave function, we start from the integral 



exp 



1 



+ a I x + fix 



H n (x)dx 



(48) 



where a is a pure imaginary quantity. As a result, we obtained [271: 



2tt 



1 - 2a 



cxp 



/3 2 



1 + 2a 



2(1 -2a) J VVl -4q 2 



Vl - 4,a 2 



(49) 



and for the evolution of the wave function 



e -i(n+l/2)t 
G(x, x't)tfi(x', t = 0)dx' = . x 

^2 n n\^pK 



x exp 



* F 2 (t) 

■ 2 ' 

sin r 



dr + x 2 (t) ctgi exp 



{x-x{t)f 



+ ip(t)x)H n (x-x(t)) (50) 



Now we use formula (44) to obtain the expression for the tomographic-probability distribution at any 



particular moment of time t from expression for the wave function ( 50 ) : 



w n (X, n, v, t) = exp 



(X - fix{t) - vp{t)f 



V 2 + fl 2 



X — pLx{t) — Vp{t) 

\Jv 2 + p 2 



(51) 



Finally, we show the second way of calculations. 



The tomographic-probability distribution at the initial moment of time t = according to formula ( 44 ) 
reads: 



w n {X,{i,v,t = 0) 



1 



1 



exp 



X 2 j Til 

,2 _i_ ,,2 ) ™ 



X 



yjv 2 + H 2 



(52) 



For the nth excited state of the harmonic oscillator, the mean values of the position and momentum are 
equal to zero, so the initial conditions for the functions x c i(t),p c i(t) are zero and these functions coincide 



with x(t),p(t). Then using formula for the evolution of the tomographic-probability distribution (40) we 



obtain the same time dependence of the tomogramm distribution (51). 



5 Optical Tomograms of the Driven Oscillator 

In this section, we present all the results obtained in the previous sections in terms of the optical 
tomogram w (X,9). The evolution equation for the optical tomogram can be derived from the Moyal 
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equation (31) using Mowing correspondence rules [28]: 

qW(q,p) i — > 



pW(q,p) i — ► 
d_ 

ai, 
a 



_d_ 

ax 

d_ 

ax 



d 



sin#— + Xcos0 w(X,t 

ao i 



-1 



cos#^ + Xsin9 | w(X,9) 

ao 



a 



. } W(q,p)^cos9— w(X,9) 
<)p W(q,p)^sme^w(X,e) 



(53) 

(54) 

(55) 
(56) 



After such substitutionswe can write Eq. (33) in the form containing the optical tomogram w(X,9,t): 

a a 



01 09 +f<t)£f**0)u<X.O.O =0 



Integrals of motion for Eq. (57) can be written in a form, similar to formulae (37)- (39): 

e (x,e,t) = e + t 

X (X, 6,t) = X - x(t) cos 9 - p{t) sin 9 



(57) 



(58) 
(59) 



In view of (58) and (59) we obtain the final expression for the time evolution of the quantum optical 
tomogram of the driven harmonic oscillator 

w(X,9,t) = wo(X -x(t) cos 9- p(t) sin 9,9 + 1) (60) 

where wq (X, 9) is optical tomogramm distribution at the initial moment of time t = and it depends on 
just two variables X, 9. 

Optical tomograms at the initial moment of time wo(X, 9) for the cases considered in previous section 

sin#. So, for the coherent state of 



cos V, v 



can be found from (|45|), (|52|) using the substitution /x 

tomogram reads 
: exp ( — {X — x c i(t) cos 9 — Pd(t) sin#) 2 



the driven oscillator the time evolution of the optical tomogram reads 

w{X,9,t) = — r 
v 7r 

In the case of the nth excited state of the harmonic oscillator, the evolution of the optical tomogram is 



(61) 



W n (X,9,t) 



1 



2 n nU/7T 



exp ( - (X - x(t) cos 9 - p{t) sin 9f ) {X - x{t) cos 9 - p(t) sin 9) 



(62) 



6 Conclusions 

To conclude, we point out the main results of our study. 

We presented a short review of the tomographic probability representation of quantum mechanics 
on the example of the driven harmonic oscillator. We considered the evolution equation for the tomo- 
graphicprobability densities of the symplectic and optical tomograms. We concentrated on the important 
solutions of the evolution equations which correspond to the Gaussian packets (normal distributions of 
the coherent state) and GaussHermite distributions corresponding to the oscillator excited states. Other 
tomographic approaches to the driven-oscillator problem like the photon-number tomographic scheme 
will be studied in future publications. 
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